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We show that scattering-based interaction among nearly degenerate optical modes is the key factor in low 
threshold generation of Kerr frequency combs in nonlinear optical resonators possessing small group velocity 
dispersion (GVD). The mode interaction is capable of producing drastic change in the local GVD, resulting 
in either a significant reduction or increase of the oscillation threshold. It is also responsible for the majority 
of observed combs in resonators characterized with large normal GVD. We present results of our numerical 
simulations as well as supporting experimental data. 
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There are several experimentally observed phenomena in 
optical frequency combs generated with monolithic microres- 
onators (Kerr combs) that have been studied by multiple re- 
searchers (see IjJ] for review), but are not as yet explained 
theoretically. An example of the observed phenomena lack- 
ing a theoretical explanation is the strong dependence on the 
selection of the externally pumped optical mode of both the 
spectral envelope and the generation threshold of the comb. 
The purpose of this contribution is to provide a theoretical 
framework for these observations. 

Generation of optical frequency (Kerr) combs results from 
modulation instability of the continuous wave (cw) light con- 
fined within a mode of the resonator 0], The comb is gener- 
ated when the power of the cw pump exceeds a certain thresh- 
old. The frequency harmonics of the comb are produced in the 
modes of the resonator adjacent to the pumped mode; these 
harmonics are phase locked 0]. 

Theoretical studies dealing with idealized nonlinear res- 
onators suggest that the comb spectra are symmetric with re- 
spect to the frequency of the optical pump HH3]. Further- 
more, the frequency spacing of the generated harmonics can 
vary from a single free spectral range (FSR) (the fundamen- 
tal comb) to multiple FSRs of the resonator, depending on the 
power and frequency of the pump as well as group velocity 
dispersion (GVD) of the resonator modes 

Theoretical predictions [4-7] describe results of experi- 
ments performed with resonators having anomalous GVD 
1 1411 . while multiple experimental observation with normal 
or nearly-zero GVD resonators JE 02] evidently contradict 
the theory. In resonators characterized with anomalous GVD 
the power-dependent frequency shift of the modes is compen- 
sated by the dispersion, leading to extended dynamic range 
for comb formation as well as ensuring excitation of the comb 
from zero field fluctuations (soft excitation regime). Theo- 
retically, there is no soft excitation regime for combs of res- 
onators with normal GVD and the dynamic range for the comb 
formation is limited there [10]. Practically, asymmetric combs 
and the soft excitation regime of the comb are frequently ob- 
served. In this work we show that the discrepancy originates 
from the interaction of resonator modes overlapping the same 
mode volume and belonging to different mode families. The 
mode interaction disturbs the regular GVD of the resonator 



spectrum so much that the local GVD significantly changes 
its value, as well as its sign. We argue that mode crossing is 
the chief reason for observation of low threshold frequency 
combs in large enough monolithic resonators. 

Very recently performed a study of universal dynamics of 
Kerr comb formation in dielectric microresonators [11]. This 
study explained the anomalously large linewidth broadening 
in octave spanning frequency combs 11211 and multiple beat- 
notes observed in low repetition rate comb systems when the 
comb was demodulated on a photodetetor 11311 . Our current 
study shows that the theoretical consideration of the low rep- 
etition rate comb systems should be adjusted to take into ac- 
count the interaction between the modes of the nonlinear res- 
onator. 

To reveal the importance of mode crossing we numerically 
simulate generation of combs in 21 identical optical modes 
coupled through cubic nonlinearity. The modes have unequal 
separation resulting from nonzero GVD of the resonator. To 
perform the simulation we numerically solve a set of nonlinear 
equations derived using the physical model presented in 12j,|7[]: 



'a 3 = -( 7o + iu)j)aj + ~[V, aj] + F e~ lult S lhj , (1) 
where V = —Tig(e' ! ) 2 e 2 /2 is the interaction Hamiltonian, 
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(aj is the annihilation operator of the field in 



the mode), g = humeri? /(Vtiq) is the interaction constant, 
u)\i is the frequency of the pumped mode, no and m are the 
linear and nonlinear refractive indexes of the material, V is the 
mode volume, 5nj is the Kronecker's delta; 70 = 70c + 7o« 
is the half width at the half maximum for the optical modes, 
assumed to be the same for the all modes involved; and 7o c 
and 7oi stand for coupling and intrinsic loss. The external op- 
tical pumping is given by Fq = (2jo c P / (Hojh)) 1 I 2 , where 
P is the value of the cw pump light. We neglect the quantum 
effects and do not take into account corresponding Langevin 
noise terms. 

We solved (T} numerically taking into account only the sec- 
ond order frequency dispersion recalculated for the frequency 
of the modes u>j as W12+W10 — 2cJn = —f32COJp SR /no, where 
^2 is the effective GVD of the modes with no interaction taken 
into account, c is the speed of light in the vacuum, u>fsr 1S the 
free spectral range of the resonator, e.g. 2ujfsr — ^12 — ^io- 
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The effective GVD value is negligible in a large WGM res- 
onator that has no mode interaction, 70 S> \D\ = |o;i2+wio — 
2w n |. For example, the mode non-equidistance can be esti- 
mated as D ~ -2?r x 39 Hz (D ~ 2tt x 1 kHz) in a 10 GHz 
FSR CaF 2 (MgF 2 ) WGM resonator pumpeded at 1550 nm. 
The value changes rather insignificantly compared with the 
typical bandwidth of the modes used in reported experiments 
(270 > 2tt x 200 kHz, or Q < 10 9 ) in smaller resonators 
(35 GHz FSR): D ~ -2vr x 5.2 kHz (D ~ 2tt x 8 kHz). 
On the other hand, D ps 270 is required for generation of the 
fundamental frequency comb with soft excitation [2]. This is 
possible in a small resonator. For instance, for a fundamental 
TE mode of a 100 GHz FSR MgF 2 WGM resonator pumped 
with 1721 nm light D « 2ir x 200 kHz. Any observation of 
soft excitation of the comb in a larger resonator is the evid- 
ance that the spacing of the modes is somehow altered com- 
pared with the value of an ideal resonator. One way for such a 
modification is related to the usage of a different mode family, 
the GVD of which can be tuned with the morphology of the 
resonator B16I1 . Another way is related to the mode interaction 
described below. 

Interaction between resonator modes has been encountered 
in WGM resonators. It is directly observed as a disruption of 
the continuity of dispersion in smaller resonators 11711 . and 
is indirectly indicated by the asymmetry of the Kerr comb 
spectra |Hl[l5t] (the spectrally narrow combs generated in res- 
onators with no mode interaction has to be symmetric because 
of energy conservation). A spectrum of a microresonator 
taken in a broad frequency range 111 811 clearly shows the pres- 
ence of spurious modes changing their position with wave- 
length with respect to the fundamental mode family. Since the 
spurious and the fundamental modes are able to interact due 
to unavoidable imperfections of the resonator shape, they also 
can change the dispersive properties of the resonator. Even if 
the spurious modes are not seen in the experiment with a par- 
ticular selection of the coupling technique, they still can exist 
in the resonator |T^, 2(1 . Only a resonator supporting a sin- 
gle mode family (see, e.g., 12111 ) can be considered free from 
mode interaction. Therefore, it is natural to expect that mode 
interaction is essential in the majority of observations of Kerr 
frequency combs in optical resonators. 

Let us assume that the mode containing the first higher fre- 
quency harmonic (012) interacts with a mode of the resonator 
(c) having the same loaded Q-factor and not belonging to the 
mode family the comb is generated in. The interaction, de- 
scribed by Hamiltonian hK(a\ 2 c + c^ayz), results in the well 
known splitting of the resonance for an (see Fig. [TJ. The 
splitting may be considered as a pure frequency shift of the 
mode a\2 expressed as 



Ul\2 — W12 
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in the asymptotic case of large difference between the eigen- 
frequencies of the interacting modes, A = uj c — W12, com- 
pared with the interaction constant k, | A| 3> k. 

The shift of mode a\i resulting from the interaction with 
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FIG. 1: Density plot illustrating frequency and amplitude depen- 
dence of the forced response of a linear mode an on the detuning 
A = u) c — W12. Mode ai2 is pumped externally with a monochro- 
matic force of a constant amplitude (|fo/7o| = 1)- The frequency 
separation between the force frequency (uj) and the frequency of the 
mode (&12) is determined as S = W12 — W12, where uj\2 corresponds 
to the frequency of the free mode. The calculations are made for 
K/70 = 20. The dependence shows that for large |A| the interaction 
primarily results in frequency shift of the mode of interest. The mode 
splits into a symmetric doublet when A approaches zero. 



mode c leads to the modification of the effective GVD value 
for the pumped mode and the first two sideband modes. Re- 
ally, using Eq. (O we find 
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According to Eq. (01, it is enough to have n 2 « — 2A70 to 
achieve the desirable value of the GVD in any resonator with 
a small intrinsic dispersion. 

To validate the analytical calculations we performed a sim- 
ulation for the 21 mode at the above conditions but also 
took into account the interaction with a mode c. Selecting 
k/7o = 20 we found that the soft excitation of the Kerr 
comb is possible for a wide range of frequency detunings: 
— O.270 > A > — 2.770 (Fig.|2]i. The Kerr comb for specifi- 
cally selected parameters of the system is shown in Fig. ([3^). 
It has a slightly asymmetric spectrum and fast roll-off of the 
higher order harmonics. The comb starts from zero fluctua- 
tions of the field with essentially zero initial conditions. Since 
such excitation regime of the comb is absent in the case where 
no mode interaction is available, we conclude that mode inter- 
action is the cause, as predicted by the reasoning above. 

According to the numerical simulations the frequency comb 
based on mode interaction has two distinct features: (i) the 
dynamic range of the soft excitation of the comb is limited 
with respect to the power of the external pump and (ii) the 
repetition rate of the comb depends on the pump power. We 
evaluated the power of the first comb sideband generated in 
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FIG. 2: Density plot showing normalized amplitude distribution of 
the first higher frequency sideband of the Kerr comb generated due 
to interaction of the resonator modes. The distribution is plotted as 
a function of the detuning between the pumping laser and the cor- 
responding mode of the resonator (Ao/70), and the frequency sepa- 
ration between the optical mode containing the comb harmonic and 
the service mode (c). Coupling parameter ft/70 = 20 is selected. 
The GVD of the resonator is selected to be D/70 = —0.02. The 
amplitude of the external cw pump is l-Fbl/70 = 2. Mode an is 
pumped. 
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FIG. 3: Spectra of the frequency combs generated in a resonator 
having nearly zero normal GVD (D/70 = —0.02) when one of the 
resonator modes is shifted due to the interaction with another mode. 
Coupling parameter K/70 = 20 is selected. The amplitude of the 
external cw pump of the central mode (an) is |Fo|/7o = 2. Three 
cases are shown: (a) mode a±2 is shifted due to the interaction; (b) 
mode ai3 is shifted; and (c) mode an is shifted. Because the relative 
GVD of the resonator is small it does not influence the nonlinear 
process, so shapes of the comb envelopes simply scale from picture 
to picture. 



the shifted resonator mode (0,12) as a function of the pump 
amplitude and frequency detuning for the fixed interaction 
value with another mode (Figj4|. The limitation occurs due 
to clamping of the value of the nonlinear frequency shift of 
the resonator modes that can be compensated by the inter- 
action of ilia and c modes. The power dependence of the 
comb frequency is illustrated by (FigO, showing the differ- 
ence [(wi2 — cjn) — (w'12 — uj)]/-fo, where u/12 is the fre- 
quency of the generated sideband. 
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FIG. 4: Normalized amplitude of the first comb sideband generated 
in the mode ai2 shifted due to the interaction with the service mode. 
Nearly zero normal GVD (.D/70 = —0.02) and coupling parameter 
K/70 = 20 are selected. 
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FIG. 5 : Frequency difference between the local FSR of the resonator 
(not modified by the interaction with the service mode) and the rep- 
etition rate of the comb as a function of the amplitude and frequency 
of the pump. 

Interaction of the modes also explains experimentally ob- 
served generation of combs having different repetition rates in 
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resonators with normal GVD fl. The comb step changes de- 
pending which of the comb-generating mode in the sequence 
is shifted due to the interaction. The intrinsic GVD of the res- 
onators is too small to influence the process. We have simu- 
lated generation of higher repetition rate combs and validated 
the concept (Fig.|3J),c). The simulation shows that the detun- 
ing of the pump light from the corresponding mode of the res- 
onator, discussed in ]8[], is of less importance for generation 
of high repetition rate frequency combs having soft excitation 
regime, compared with mode interaction considered here. The 
change of the comb step in normal GVD resonators is differ- 
ent from the mechanism of the step change in resonators with 
anomalous GVD, where repetition rate strongly depends on 
the intracavity power of the pump light |0, [3 . 

We performed an experiment with a large WGM resonator 
having nearly zero relative GVD (D/70) and observed gen- 
eration of frequency combs with envelopes similar to those 
predicted by the theory. In our experiment we used a CaF2 
resonator having 6721 fim in diameter. The resonator had 
approximately 9.9 GHz FSR with loaded quality factor ex- 
ceeding 10 9 . We pumped the resonator with 1545.5 nm light 
emitted by a distributed feedback semiconductor laser. The 
light was coupled to the resonator via a coupling prism. The 
optical power emitted by the laser was 15 mW, and 3.2-1 mW 
of the light entered the selected modes of the resonator (the 
value depends on the selected mode). The output light was 
collected using a PM Panda fiber and introduced to an optical 
spectrum analyzer. The resultant optical spectra are shown in 
Fig. dill- 
Frequency combs simulated numerically have frequency 
shapes similar to the combs observed experimentally (com- 
pare Fig.[3]and|6]l. The comb properties change significantly 
when we pump different modes of the resonators. Such a 
modification cannot be explained by the change in the geo- 
metrical part of the GVD of the modes. The interaction with 
the degenerate modes, on the other hand, perfectly explains 
the observation. 

To conclude, we have shown theoretically that experimen- 
tally observed generation of optical Kerr frequency combs in 
nonlinear resonators with small group velocity dispersion re- 
sults from the linear interaction of resonator modes. The mode 
interaction changes the frequency of the modes enabling soft 
excitation regime of the frequency combs. Only combs pro- 
duced via hard excitation can otherwise be generated in these 
resonators. 

The authors acknowledge partial support of the reported 
study by DARPA's IMPACT program. 



[1] T. J. Kippenberg, R. Holzwarth, and S. A. Diddams, Science 
332, 555 (2011). 

[2] A. B. Matsko, A. A. Savchenkov, D. Strekalov, V. S. Ilchenko, 
and L. Maleki, Phys. Rev. A 71, 033804 (2005). 

[3] P. Del'Haye, A. Schliesser, O. Arcizet, T. Wilken, R. 
Holzwarth, and T. J. Kippenberg, Nature (London) 450, 1214 




PL 1 r »M . 1 , Tw™ f« 30 

1547.5 1548.0 1548.5 1549.0 
Wavelength, nm 



FIG. 6: Experimental spectra of three optical frequency combs gen- 
erated in the same overmoded WGM resonator when light is cou- 
pled to three arbitrary modes having strong interaction with other 
modes. The wavelength of the pumping light was changed to se- 
lect the modes, and depending on the mode selection we observed 
(a) single-FSR, (b) dual-FSR comb, and (c) triple-FSR combs. Gen- 
eration of higher order frequency combs, not shown here, was also 
recorded. 



(2007). 

[4] I. H. Agha, Y. Okawachi, M. A. Foster, J. E. Sharping, and A. 

L. Gaeta, Phys. Rev. A 76, 043837 (2007). 
[5] I. H. Agha, Y. Okawachi, and A. L. Gaeta, Opt. Express 17, 

16209 (2009). 

[6] Y. K. Chembo, D. V. Strekalov, and N. Yu, Phys. Rev. Lett. 104, 
103902 (2010). 

[7] Y. K. Chembo and N. Yu, Phys. Rev. A 82, 033801 (2010). 

[8] A. A. Savchenkov, A. B. Matsko, V. S. Ilchenko, I. Solomatine, 
D. Seidel, and L. Maleki, Phys. Rev. Lett. 101, 093902 (2008). 

[9] O. Arcizet, A. Schliesser, P. DelHaye, R. Holzwarth, and T. J. 
Kippenberg, "Optical frequency comb generation in monolithic 
microresonators," in Practical Applications of Microresonators 
in Optics and Photonics, A. B. Matsko, ed. (CRC Press, 2009), 
Chap. 11. 

[10] A. Matsko, A. Savchenkov, W. Liang, V. Ilchenko, D. Seidel, 
and L. Maleki, "Group Velocity Dispersion and Stability of 
Resonant Hyper-Parametric Oscillations," in Nonlinear Optics: 
Materials, Fundamentals and Applications, OSA Technical Di- 
gest (CD) (Optical Society of America, 201 1), paper NWD2. 

[11] T. Herr, J. Riemensberger, C. Wang, K. Hartinger, E. Gavartin, 
R. Holzwarth, M. L. Gorodetsky, and T. J. Kippenberg, 
\arXiv : 1111.3071i ;l. 

[12] P. Del'Haye, T. Herr, E. Gavartin, M. L. Gorodetsky, R. 
Holzwarth, and T. J. Kippenberg, Physical Review Letters 107, 
063901 (2011). 

[13] S. B. Papp, S. A. Diddams, Phys. Rev. A 84, 053833 (2011). 
[14] W. Liang, A. A. Savchenkov, A. B. Matsko, V. S. Ilchenko, D. 

Seidel, and L. Maleki, Opt. Lett. 36, 2290 (201 1). 
[15] I. S. Grudinin, N. Yu, and L. Maleki, Opt. Lett. 34, 878 (2009). 
[16] A. A. Savchenkov, A. B. Matsko, W. Liang, V. S. Ilchenko, D. 

Seidel, and L. Maleki, Nature Photonics 5, 293 (201 1). 



5 



[17] P. Del'Haye, O. Arcizet, M.L. Gorodetsky, R. Holzwarth, T.J. 

Kippenberg, Nature Photonics 3, 529 (2009). 
[18] F. Ferdous, H. Miao, D. E. Leaird, K. Srinivasan, J. Wang, L. 

Chen, L. Tom Varghese, and A. M. Weiner, Nature Photonics 5, 

770(2011). 

[19] A. A. Savchenkov, A. B. Matsko, V. S. Ilchenko, D. Strekalov, 
and L. Maleki, Phys. Rev. A 76, 023816 (2007). 



[20] T. Carmon, H. G. Schwefel, L. Yang, M. Oxborrow, A. D. 

Stone, and K. J. Vahala, Phys. Rev. Lett. 100, 103905 (2008). 
[21] A. A. Savchenkov, I. S. Grudinin, A. B. Matsko, D. Strekalov, 

M. Mohageg, V. S. Ilchenko, and L. Maleki, Opt. Lett. 31, 1313 

(2006). 



